The low-T phase diagram of LiHoxYi_xF4 
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The LiHoxYi_xF4 compound is widely considered to be the archetypal dipolar Quantum Ising 
system, with longitudinal dipolar interactions Vif' between Ho spins {i, j} competing with transverse 
field-induced tunneling, to give a T = quantum phase transition. By varying the Ho concentration 
X, the typical strength Vb of V^" can be varied over many orders of magnitude; and so can the 
transverse field H±^ . A new effective Hamiltonian is derived, starting from the electronuclear degrees 
of freedom, and valid at low and intermediate temperatures. For any such dipolar Quantum Ising 
system, the hyperfine interaction will dominate the physics at low temperatures, even if its strength 
Aq < Vq: one must therefore go beyond an electronic transverse field Quantum Ising model. We 
derive the full phase diagram of this system, including all nuclear levels, as a function of transverse 
field H±, temperature T, and dipole concentration x. For LiHoxYi_xF4 we predict a re-entrant 
critical field as a function of x. We also predict the phase diagram for x= 0.045, and the behavior 
of the system in magnetic resonance and /^SR experiments. 
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I. INTRODUCTION 

A. The TFQI Model for LiHoxYi-xF* 

For at least a decade the LiHoxYi_xF4 compound has 
been considered to be an ideal experimental reahzation 
of the well-known 3-dimcnsional transverse field Quan- 
tum Ising model (TFQIM). According to this view, at 
temperatures well below an anisotropy energy fio, it is 
described by the Hamiltonian: 



(1) 



where Tj is a Pauli vector describing a two-level effec- 
tive electronic spin at spatial position r = r^, Vf^'^ is a 
longitudinal inter-spin interaction, with nearest-neighbor 
strength J7o which, depending on the dilution x, can have 
either ferromagnetic or a frustrating character; and the 
'transverse field' term Aq '''i is controllable externally 
(usually by applying a transverse magnetic field). The 
most distinctive feature of the TFQI Model ([T|) , which is 
central to the whole field, is the competition between Vb, 
which tries to order the system, and Aq, which causes 
quantum fluctuations out of the ordered state. At T = 
one expects a quantum phase transition between ordered 
and quantum disordered states when Aq/Vq ^ 1, and this 
is probably the simplest theoretical example of a quan- 
tum phase transition. The apparent confirmation of this 
'quantum critical' picture for LiHoxYi_xF4 has lent con- 
siderable importance to the experiments on this system. 

The main arguments in favor of this picture for 
LiHoxYi_xF4 are as follows: 

(i) The strong crystal field Ho single-ion anisotropy 
yields an Ising doublet ground state, with a crystal field 
Hamiltonian yielding an appreciable Aq at small H^. 
The dominant inter-Ho spin-spin interaction is dipolar, 
with strength Vb(x) = o;^; with a ^ 1 in 

Kelvin units. Thus when x= 1 one expects a dipolar- 
ordered ferromagnetic phase below ^ IK, which is ob- 



served; it exhibits both classical and quantum phase tran- 
sitions to the paramagnetic phase [ij. 

(ii) In LiHoxYi_xF4 the magnetic Ho ions and the non- 
magnetic Y ions have very similar atomic volumes; dilu- 
tion of the Ho by Y is then possible with negligible distor- 
tion of the lattice. This dilution weakens the interactions 
and introduces randomness and frustration; very differ- 
ent physical regimes can then be studiedjl, [1]; see Fig. 
1 in Ref. 0- In particular, one expects a low-T spin- 
glass phase at small x, below a transition temperature 
Tc ^ ax. At x=: 0.167 a spin-glass phase is foundjl, Q at 
low T and H±^ — 0, with a crossover to the paramagnetic 
phase at higher T and H±. At x= 0.44 the tunneling 
of domain walls in the ferromagnetic phase was found [6| 
and differences between quantum and classical anneal- 
ing protocols were observed^- At x= 0.045 the system 
shows a peculiar narrowing of the spin fluctuation spec- 
tral width as temperature is decreased described as 
"anti spin-glass" behavior. 

(iii) For extreme dilution one expects single Ho ion 
behavior. In experiments at x= 0.002, hysteresis loops 
of the magnetization due to single spin tunneling are 
observed [9j (co-tunneling of pairs of spins was also ob- 
served at x= 0.002, showing that interaction effects can- 
not be neglected even at this dilution [13, [HI)- 

Thus, according to these arguments, a TFQI model 
like ([T|) should describe LiHoxYi_xF4 for all x, provided 
kT,ijLBH± ^ r^o; and as such, LiHoxYi_xF4 should be 
a model system for all dipolar magnets. However we 
argue in this paper that the LiHoxYi_xF4 system (and 
by implication, many other dipolar magnets) need to be 
described in a quite different way. There are two main 
problems with the simple TFQI picture, both noted and 
analysed in ref.[l3|- These are 

(a) Hyperfine Interactions: The on-site Ho hyperfine 
interaction is not small - in fact even at x— 1, 
Aq ^ Vb, and for x^ 1, the hyperfine interaction is over- 
whelmingly dominant! A few experimental papers have 
heeded this point, remarking (i) that even the x= 1 phase 
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diagram, near the T = ferromagnetic-paramagnetic 
transition, is modified by the hyperfine interaction[l|; 
and (it) that the nuclear spin bath, considered now as a 
quantum environment [isj . should strongly affect the Ho 
spin dynamics near this quantum critical point p^. ITU- 
However we shall show here that the effect of nuclear 
spins on dipolar magnets is much more profound than 
this, even when the hyperfine interaction is quite weak. 
This very surprising result means that one must recon- 
sider the application of the TFQI Hamiltonian to a large 
variety of systems, hitherto analysed without reference 
to the hyperfine couplings. 

(b) Transverse Dipolar Interactions: When 1 these 
interactions add a quite large contribution to the trans- 
verse field - to quantitatively understand the phase di- 
agram one then needs to include themfl^, [l^ (see also 
ref.p7l|), both in the spin-glass and in the ferromagnetic 
regimes 111 M M MM ■ 



B. An ENQI Model for Dipolar Ising Magnets 

To properly treat the physics of Quantum Ising sys- 
tems, we have to recognize that the use of a simple pa- 
rameter Ao(i?_L); introduced a long time ago by exper- 
imentalists as a convenient way of defining an effective 
transverse field acting on the Ising spins, is actually mis- 
leading. Because of the nuclear spins, the true effective 
transverse field in a quantum Ising system is very differ- 
ent from Aq; moreover it depends on the actual nuclear 
spin state of the system. 

In what follows we will derive a theoretical framework 
with the nuclear spins included from the beginning. The 
system is described at low energies in terms of 'electronu- 
clear' complexes which interact via renormalized dipolar 
interactions. In its general form [see Eq.® just below] 
this 'Electronuclear Quantum Ising" (ENQI) Hamilto- 
nian includes all the nuclear spin levels. However at very 
low T or for small x, we can use a much simpler Hamilto- 
nian referring only to the lowest electronuclear doublet, 
and this takes the form 

H=-Y.Vi^\H^)sls]~K{H^)Y,^'i. (2) 

where now Sj operates only on the single electronuclear 
doublet involving the nuclear states with = ±/. Now 
this simplified model looks like the standard TFQI model 
in (IT|), but it behaves very differently - both Y{!j''{H±) 

and A(i7j^) are renormalized from their original values 
in ([T]), and they depend strongly on H±^ [in the case of 
A(i7j^), this dependence is radically different from that 
in the original parameter /S.q[H±_)\. The strength and 
behavior with field of these variations depends crucially 
on the strength of the hyperfine interaction; moreover, 
as noted above, we must use this ENQI model at low T 
even when the hyperfine coupling Aq -^Vq^ which is more 
typical for a general anisotropic magnet. 



More generally, when kT is not small compared to 
the splitting between nuclear levels, we must define a 
set of 2/ + 1 electronuclear 'pseudospins' (each of which 
are spin- 1/2 doublets) labelled by quantum numbers 
m = 1,1 — 1, ... — /, an occupation number riim for the 
occupation of a given pseudospin on site i, and a set of 
pseudospin operators Sim and pseudospin energies e™; we 
have the general ENQI Hamiltonian 

- 'Et,,n nrm[e,n + A,n{H^)sfJ , (3) 

where the V^jfmm'(-ff±) represent interactions between 
pseudospins m,m' on different sites and the transi- 
tion matrices A^ only operate on individual pseudospins, 
ie., within the space of each electronuclear doublet on a 
given site. We can think of a set of 2/ + 1 independent 
quantum Ising systems, each having a different 'trans- 
verse field' Am, which however can interact via the lon- 
gitudinal fields Vfj^j^^, (iJ j_ ) . 

In disordered dipolar-coupled spin systems, one must 
also add a term which describes the random transverse 
couplings in the system. Its detailed form is given in 
section IIIIl and its quantitative effects are discussed in 
section IIVI 

In this paper we concentrate on the LiHoxYi_xF4 sys- 
tem, for which precise results and experimental predic- 
tions can be established for the phase diagram, so it 
can be used as a test case. The effective Hamiltonian is 
strictly applicable to systems where Aq,Vo fioi in the 
regime where T,^^H±^ <^ fig- This approach enables (i) 
illumination of the relevant physics of the LiHoxYi_xF4 
system (ii) generalization to other systems, e.g. sys- 
tems in which <C Vq (see Sec. IIVP (in) construct- 
ing a framework for the treatment of dynamical prop- 
erties. However, in the LiHoxYi_xF4 system the con- 
dition AojVo ^ i^o is not that well satisfied, and the 
condition T,hbH± <^ flo, while is satisfied in the whole 
relevant phase diagram at low x, is not satisfied near 
criticality at large concentrations. For this reason, and 
since single ion properties dictate much of the physics in 
the LiHoxYi_xF4 system, we also use exact diagonaliza- 
tion of the Ho electronuclear spin states. This enables us 
to give quantitative predictions regarding the single ion 
characteristics, and with the use of mean field approxi- 
mation, to predict the form of the phase diagram for all 

X. 

Most of the results here are new. We analyze in de- 
tail the form of the electro-nuclear states of the single 
Ho ion as function of H±, and its consequences in terms 
of entanglement entropy and magnetic resonance exper- 
iments. We show that the peculiar crystal field Hamil- 
tonian of LiHoxYi_xF4 results in a well defined Ising 
system even at high transverse fields (where Ising sym- 
metry is usually destroyed) . We obtain a general effective 
Hamiltonian valid for thermodynamic properties, incor- 
porating all 16 low e nerg y states, therefore generalizing 
the treatment in Ref.[l3] to the regime Aq < T <$:. flo- 
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We give a discussion of the phase diagram for general 
concentration x, temperature T, transverse field H± and 
hyperfine coupling Aq. With relevance to general mag- 
netic systems we show that the hyperfine interactions 
dominate the physics at low T even when Aq <^ Vq. By 
comparing the phase diagrams at x= 0.045 and x= 0.167 
we predict a novel reentrancc of the crossover transverse 
field between the quasi-SG and PM phases at low T as 
function of x, resulting from the interplay between the 
hyperfine and off-diagonal dipolar interactions. We then 
add a novel perspective to the unsolved question of the 
nature of the low-T phase at x= 0.045. Finally, we dis- 
cuss further experimental consequences of our results. 

The paper is organized as follows: In Sec.|TT]the various 
terms in the microscopic Hamiltonian for LiHoxYi_xF4 
are introduced and quantified, and single ion properties 
are analyzed. In Sec. IIIII the full low energy effective 
Hamiltonian is derived, including the transverse hyper- 
fine interactions and the off-diagonal terms of the dipolar 
interaction. In Sec. IIVI we obtain the phase diagram of 
the LiIl0xYi_xF4 system at different dilutions; we obtain 
quantitative agreement with the experimental phase dia- 
gram at x= 0.167, make predictions regarding the phase 
diagram at x= 0.045, and discuss the nature of the low 
temperature phase. In Sec. |V] we suggest experiments 
that can directly check our theory, and in Sec. IVII we 
state our conclusions. Some details regarding the deriva- 
tion of the effective Hamiltonian and the calculation of 
the phase diagram in mean field are deferred to appen- 
dices. 



II. INTERACTIONS IN THE LiHoxYi_xF4 
SYSTEM 



In this section we give the quantitative form of the 
LiHoxYi_xF4 Hamiltonian, which is a sum of crystal 
field [qI. [2H. Zeeman, inter-Ho, and hyperfine interaction 
terms: 



A. TFQIM terms 

Let us first consider the terms which feed directly into 
the TFQIM Hamiltonian ([1]), ie., the terms Hcf, Hz, and 
Hint- For LiHoxYi_xF4 these are given in turn by: 

(i) The 'crystal field' term Hd includes the single-ion 
crystal field and spin-orbit terms [2l|. Because of the 
very strong spin-orbit coupling, J is a good quantum 
number for the Ho ion with J = 8. A crystal field term 
of form (J^ + Jt) strongly mixes states with differing 
by ±4 [3, H^, and a strong term severely distorts 
the level spacing. There are other terms as well - for 
computations in this paper we will use a form written in 
terms of the usual Stevens operators as[l^. [2l]| 



H 



cf 



i=2,4,6 



B?0^ 



BlOtiS) 



1=4,6 



BfofiC) (5) 



with values assumed to be l2l| 



= -0.696, b2 = 4.06 x 10 



,-3 



= 4.64 X 10"^, Bi{C) = 4.18 x 10-^ 
B^{C) = 8.12 X 10^"^, B^{S) = 1.137 x 10"^ 



(6) 



from which we sec that the Oq{C) term also has a non- 
trivial effect, bringing in a (J^ + Jl) term (note that to 
properly judge the relative importance of coefficients _B™ 
and B™ with / 7^ I', we should directly compare ,PB™ 
and J'-' B\]\ and not B™ and B™). 

The ground state is an Ising doublet, with states de- 
noted here by W) and ||), which mix states with Jz = 
±7, ±3, ^1,^5. The first excited state | is roughly 
rio = 10.5i4r above the ground state doublet, and is a 
mixture of Jz = 6,2,— 2,— 6. The other 14 states are 
much higher in energy, and the total span of the J = 8 
manifold is roughly df = 500if [Hi. 

(ii) The Zeeman coupling to the Ho spins is given by 
the usual form 



Hz = - ^gjUBH ■ Ji 



(7) 



H — Hcf + Hz 



H 



int 



H 



hyp ■ 



(4) 



Note that we have dropped: (i) the spin-phonon interac- 
tion, important for spin relaxation (ol. [lOl. [2^ : (ii) hyper- 
fine interactions between the Ho ion and other nuclear 
species (F, Li) as well as with Ho nuclei on nearby sites; 
and (iii) the nuclear Zeeman couplings. None of these 
terms have an appreciable effect on the phase diagram of 
LiHoxYi_xF4. Note however that they will be crucial for 
the low-T Ho spin dynamics, since even very small hy- 
perfine terms can strongly affect relaxation dynamics and 
decoherence in the low-T quantum regime[23l. [23|. where 
phonon relaxation is also important in strong transverse 
fields [2i,[2i. 



with 9,7 = 5/4. We are particularly interested in the ef- 
fect of a transverse field H± <C fio/A'B, which induces a 
coupling Aq between the two Ising ground states in sec- 
ond order perturbation theory via the state | Fj). Thus, 
for small fields Aq oc i?^ ; by putting in the numbers one 
finds 



Ao(i/±) ~ 9(mbB±)V^^o « 0.4T2[K] 



(8) 



in Kelvin units (see e.g. Figs. 1,2 in Ref. [2l|). At larger 
fields, H± « f2o/^B ~ 2T, perturbation theory breaks 
down, I Fj) mixes strongly with |t) and |i)[2l|, and Aq 
is approximately linear in H±. An important feature of 
LiHoxYi_xF4 is that ^ ^o- Thus, the system stays 
Ising-like even when H± > Qo/fiQ, deep inside the para- 
magnetic regime (see FiglT|). This contrasts with most 
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FIG. 1: Energy of the 17 electronic states of a single Ho 
ion in LiHoxYi_xF4, ignoring all inter-Ho interactions, as a 
function of transverse field H±. The Ising-like character is 
well maintained until H± ~ 20T. 



other anisotropic dipolar systems, which are dominated 
by easy- axes terms, so that the same energy scale dic- 
tates the anisotropy and the quantum fluctuations, and 
in the quantum phase transition regime H± w there 
is no real Ising character. 

(iii) The intcr-Ho ion spin-spin interactions depend 
strongly on the Ho concentration x. They have the gen- 
eral form 



Hi, 



(9) 



Experiments [27| and theoretical analysis [21| both show 



that U^j^ is dominated by the dipolar interaction, ic. 



a/3 
ij 



K 



a/3 



a/3 

ij ~ J2 ' 

|r»i|5 



(10a) 
(10b) 



Here the strength of the nearest-neighbor dipole-dipole 
interactions between the spins 3i , 3j is 



Un 



47r Vr. 



(11) 



where Vc is the unit cell volume; the unit cell size is 
(1,1,2.077) in units of a = 5.175A, with, when x= 1, 
four Ho ions per unit cell, at positions (0,0,0), (0, a/2, 
c/4), (a/2, a/2, -c/2) and (a/2, 0, -c/4); and r.y n-rj. 
When x< 1 some of the Ho ions are substituted by Y ions, 
and the couplings acquire a random distribution whose 
character depends strongly on x. The typical value Uq 
of the nearest-neighbor coupling then becomes roughly 
IJo ~ xUo- For LiHoxYi_xF4, Uq ^ 0.3K. Note however 



that the energy Vq characterizing the total effect of the 
longitudinal dipolar interactions is somewhat larger than 
Uq, since as noted in the introduction, Vo(x) = J2j{^iT'^- 
Thus the strength of Vq depends on how the spins are 
arranged. It can be estimated from the ordering temper- 
ature, and typically Vo/Uq ^3 — 5. One can see depar- 
tures from linearity in x; for example in the diluted sys- 
tem, even for rather small x, close pairs and even triplets 
can dominate certain properties. There are also antifer- 
romagnetic exchange interactions between the Ho ions, 
which for x= 1 were measured to be about half of the 
nearest neighbor dipolar interaction 27| . Therefore, the 



exchange interactions have little quantitative significance 
even for the undiluted LiHoFzt |21| . and are completely 
negligible for x^ 1. 

If we now take these 3 terms and truncate the Ho ions 
to their lowest doublet, we get back the TFQIM in ([1]), 
which predicts a quantum phase transition for x= 1 at a 
transverse field where Aq ^ Vq, ie., at H± 3T[iI. \2l\. 
In fact the actual transition happens at i?^(x = 1) = 
4.9T[i1, [2l| , which is the first sign that there is something 
wrong with this naive picture. To see what is going on 
we now have to include the hyperfine coupling. 



B. Hyperfine interactions 

The hyperfine coupling of a single Ho atom with its 
own / = 7/2 nuclear spin gives a term 



hyp 



Jr 



(12) 



with Aj = O.O39K0. Here we ignore quadrupolar terms 
as well as the hyperfine interactions to all other species 
(Li,F, and other Ho ions); both are an order of magnitude 
smaller [28j. and hardly influence the phase diagram. 

At low energies, in the lowest doublet states |T),|i)j 
the longitudinal hyperfine term Hf^yp = AjPJ^ splits 
each electronic state into an eightfold multiplet of nearly 
equidistant levels, with separation ~ 205mi^[l| between 
adjacent levels, ie., we can write 



TTZZ 



(13) 



where f operates on the electronic doublet and for 
LiHoxYi_xF4, ujQ ^ 205 mK; this corresponds to a spin 
moment (J^) ~ bfiB for the lowest doublet. 

One can see without any reference to experiments on 
the phase diagram that the TFQI model cannot possibly 
be right at low transverse fields, using Fig[21 The lowest 
energy Ising doublet states a, a have a definite nuclear 
spin (/^ = -7/2 for |t), and 4 = 7/2 for when Hj_ = 
0. A transverse magnetic field couples a = |t,— 7/2) to 
b = li, -7/2) and d = |i, 7/2) tob= |t, 7/2), and cannot 
induce quantum fluctuations between the relevant Ising 
doublet ground states at all, but only renormalizc their 
effective spin. Only the transverse hyperfine term ^^hyp = 



5 



b b 




I) by the longitudinal hyperfine interaction. The ground state ^ 



doublet states a and d have a definite and opposite nuclear 
spin, ±7/2. A transverse magnetic field H± couples states 
with the same nuclear spin, as shown by the dashed lines. 



Aj{I^J +1 J+)/2 can change I^, and allow transitions 
between the Ising doublet states; but this hardly operates 
if ijbH^ < flo- 

Thus hyperfine interactions must be included in any 
truncation of the system to a low-energy Hamiltonian. 
Their general effect is to suppress quantum effects at low 
fields. We shall see that they are important even when 
Aq <C Vq (note that the simple argument above, showing 
the importance of the hyperfine effects, makes no refer- 
ence to the strength of these interactions!). 

C. Single Ho ion - exact results for low energies 



-2 



2 4 

(Tesia) 

FIG. 3: The 16 lowest electro-nuclear energy levels of 
LiHoxYi-xF4, plotted as a function of transverse magnetic 
field H±. The zero of energy is defined by the (field depen- 
dent) mean of levels 8 and 9. Symmetric/antisymmetric states 
are plotted in solid green/dashed blue lines. 



For H± / {p-b{J z)) , H^y^ mixes appreciably elec- 

tronuclear states with different values of Iz- This is best 
seen by performing an exact diagonalization of the full 
single Ho Hamiltonian H = Hd + Hz + -ffhyp in the 136 
eigenfunction space (17 crystal field x8 nuclear states). 
In Fig. [3] we plot the spectrum of the lowest 16 levels, 
corresponding to the electronic ground state doublet, as 
function of H±. Most generally, each of the 16 states 
can be written in the form ^oiMm \ M,m), where 
M{m) denote the z component of the electronic (nuclear) 
spin. Plotted in solid line are symmetric eigenstates, 
with aMm = oi-M-im and in dashed line anti-symmetric 
eigenstates, with aMm = —ot-M-m- 

At low fields, the electro- nuclear entanglement is 
strong, and states are given, to a good approximation, 
by the form in Ea. (|A3|) . One can then define the split- 
ting between each pair of time-reversed states by A,„, 
which are plotted in Fig[4l We findp^. as we expect, 
that A7/2 is small up to Hi^ « 2T, at which point 
(t|/^B^f±| Fj) ~ 51o- Am increases more rapidly as \m\ 
decreases, simply because for smaller |m|, transitions 
between the 2 low-energy time-reversed states can be 
achieved by lower orders in perturbation theory in H^^p, 
one then sees appreciable coupling at lower H±. As H± 
continues to increase. Am increases rapidly and eventu- 
ally saturates at a field H*^{m). Note however that the 
spectrum in Fig[3]is not symmetric. This is because tun- 
neling between the lower pairs is allowed via the state 
I F2) at energy r^o- whereas tunneling between the upper 
pairs must involve the higher excited states, at energy 



EcF higher than the lowest states. Consider, eg., the 
pairs IT, -1/2), li, 1/2) and |T, 1/2), |i, -1/2). The first 
pair has a finite matrix element in second order perturba- 
tion (T,-l/2 I if. J.I F^)(F^ |/,J, |i,l/2), which is first 
order in iij_, and gives a low-field splitting (x H^Ix/^Oi 
the second pair has a term of similar form which however 
passes via the states in the multiplet Jz = 8, 4, 0, —4, —8, 
and so gives a low-field splitting cx HxIx/Ecf, roughly 
an order of magnitude smaller. 

For H± > H]^{m) different values of m are well 
mixed, the electron and nuclear spins get disentangled, 
and the spectrum separates to two groups of 8. For 
Hx_ > H*[{7/2) the eigenstates can be approximated by 
I "0/) I V",/)- The electronic state hybridizes strongly the 
level I F2) with the ground state doublet. The states 
in the bottom group are approximately symmetric with 
respect to the electronic degrees of freedom, i.e. have 
ctMrn ~ ct-Mrn, whilc the statcs in the upper group 
have ~ ~ct~Mm- In each group states separate 

to pairs of symmetric and antisymmetric states, as noted 
above. For large H± the lower level of each pair has 
ctMm ~ dM-m and the higher level has aMm ~ —ctM-m- 
Both the energy spectrum and the form of the eigenstates 
discussed above should be revealed in electromagnetic 
resonance experiments. In Sec. IV AH we give predictions 
for such possible experiments, and their relation to the 
calculated entanglement entropy. 
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(Tesla) 

FIG. 4; The quantum fluctuation amplitudes Am induced by 
H±, plotted for \m\ = 7/2, 5/2, 3/2, 1/2; the dashed line shows 
the effective dipolar interaction V^^"^ at x= 0.167 for |m| = 
7/2 (the low-T phase transition occurs when this interaction 
~ A7/2). Note how small is A„i for large |m| and small H±. 



III. LOW-T ENQI EFFECTIVE HAMILTONIAN 

We now incorporate all the terms in with all nu- 
clear levels and the off-diagonal dipolar interactions, into 
the full ENQI model for the LiHoxYi_xF4 system, in- 
cluding all terms relevant to the phase diagram at ener- 
gies > lOmK. 

We begin by dividing up the original Hamiltonian (j4|) 
into the form 



HI" 



H 



(14) 



where 



TJZZ 

hyp 



Hz 



H"- = U, 



H^ = 



dip 



(15) 



and where we have written the dipolar interaction in the 
form 



U. 



Q/3 



(16) 



with a non-diagonal term 



^7^- 

J2 V 
J2 ^'^ij 



where TZff was defined in (fTOb]) 



(17) 



In App|X] we derive a low energy effective Hamilto- 
nian valid for T ^ fig, Mb^j_ ^ ^^o- We do this in 
three steps. We first derive the effective Hamiltonian 
for H = Hq + HI", including Ising interactions terms 



We then add -ff^yp which 



only, and obtain Ea. (jA13p 
introduces a quantum term and obtain Ea. (|A22p . We fi- 
nally include U^^p, which introduces an effective random 



field 16 1 , and an enhancement of the effective transverse 
field 1131, and as a final low-T effective Hamiltonian for 



the LiHoxYi_xF4 system we obtain 



H, 



E 



V"" ,(H^ H-^] n n- / , 



+ A,„(i// 



(18) 



Here A,„ arc the effective transverse fields acting on time 
reversed states with a given |m|, as defined in Eq. (|A10p 
(see FigU), and 7f is an effective random field, defined 
in Eq. (|A29p . We note expHcitly the dependence of the 
interactions and the effective transverse fields and ran- 
dom field on the site-dependent total transverse field H,]^ 
(IM) . 

In the low-T limit kT < coq, and for iJj^ < Hl{7/2), 
we obtain 



- Y,K{Ht)st + Y.^t{H±)st. 



(19) 



This Hamiltonian applies for any x< 1, irrespective of 
what thermodynamic phase results from it. Thus in 
LiHoxYi_xF4 it is valid for both the SG and FM regimes 
(the X dependence enters in the interaction terms and in 
the effective fields). Note, that A and 7f have a very dif- 
ferent dependence on H± and dilution x. Thus, in the FM 
phase A and 7^ are independently tunable, by changing 
X and H±_ In LiHoxYi_xF4 one may thereby real- 

ize, for the first time, both the quantum and the classical 
random field Ising models in a ferromagnetic system (see 
the theoretical details in Ref . [20| . and the experimental 
realization in Ref. [2^). 

We emphasize the essential role played here by the nu- 
clear spins. They block quantum fiuctuations. This is 
especially important for the LiHoxYi„xF4 system, whose 
peculiar crystal field Hamiltonian allows electronic tun- 
neling at second order in H±. If we drop the nuclear 
spins we can get erroneous results (eg., that the effec- 
tive random field must come at the expense of appre- 
ciable quantum fluctuations [Tsj ) . For some purposes 
one can circumvent a proper treatment of the hyper- 
fine interactions by considering a simplified crystal field 
Hamiltonian p^. [l9l. [20j (see also [lOl); where tunneling 
between the electronic spins is in high order perturba- 
tion. This gives the correct effective random field, and 
the re-entrance of the cross-over H±^ as a function of dilu- 
tion (see Sec lIVC|) . However, for other purposes a proper 
treatment of the hyperfine interactions is essential - eg., 
for the temperature and field dependence of the phase 
diagram (see Sec lIVp . and for all of the dynamic proper- 
ties. 
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IV. MEAN FIELD TREATMENT OF THE 
PHASE DIAGRAM 



The phase diagram of Quantum Ising systems Uke 
LiHoxYi_xF4 has been the object of extensive study for 
over three decades, and it was reahzcd early on that 
strong hyperfine interactions might be important [3l|. In 
the case of LiHoxYi_xF4, for x= 1 the phase diagram 
was calculated in mean- field [H, including the longitudi- 
nal hyperfine interactions; this gave an enhancement of 
the critical transverse field at low temperatures. 

In this section we analyze the phase diagram of the 
LiHoxYi_xF4 system in various dilutions, where disor- 
der effects have to be accounted for. We first discuss the 
phase diagram of a model Hamiltonian, including first the 
longitudinal dipolar and hyperfine interactions, and then 
adding the transverse hyperfine interaction. Using this 
model we make predictions for the behavior of the phase 
diagram and of the magnetization of a general anisotropic 
dipolar system where the conditions Aq,Vo <C flo are 
well satisfied, for an arbitrary ratio of Aq/Vq. This anal- 
ysis also pinpoints the basic physics dictating the phase 
diagram in the LiIl0xYi_xF4 system. However, in the 
LiHoxYi_xF4 system the condition Aq <C flo is not that 
well maintained. To solve for the phase diagram of the 
LiIl0xYi_xF4 system we use exact diagonalization of the 
single Ho ion, a mean field approximation for the inter-Ho 
interactions, and we take into account the enhancement 
of the effective transverse field by the transverse dipolar 
terms. We shall see that this then gives very accurate re- 
sults for the phase lines for LiHoxYi_xF4 when x= 0.167, 
and we make predictions for x= 0.045. 

Finally, we discuss the nature of the phases at low T; 
this is currently rather controversial. The hyperfine in- 
teractions again play a central role, in reducing quan- 
tum fluctuations and slowing the relaxation of the sys- 
tem to equilibrium in the low-T quantum regime. At 
finite transverse field we discuss the effect of the effec- 
tive longitudinal random field, emerging from the applied 
transverse field. 
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FIG. 5: Plot of the phase line separating the ordered and 
the paramagnetic phases for the model (|A2|I . in the regime 
where Aq 3> Vo, assuming Vo = l,Ao — 35, and Qq ~ 100 
in arbitrary units. The dashed line shows Eq. (|20|) . plotted 
neglecting transverse hyperfine interactions; we see that Ac 
diverges as T — > 0. The blue solid line describes the deviation 
from the classical Ising model at high transverse fields, where 
the transverse hyperfine interaction are significant, giving a 
QPT at Ao ~ flo (see Sec lTVB| . 



1. Strong hyperfine interactions 

When Aq > Vq, and kT < ajQ, the relevant Hilbert 
space comprises the lowest two electronuclear Ising-like 
levels, and we consider the Hamiltonian (|A2p . which re- 
duces to the classical Ising Hamiltonian ijjjj given in 
(jA9[) . The only effect of H± is to renormalize Fj^^ to 

Vfj^ = rf'V^Y ■ We can then immediately deduce the 
whole phase diagram of the system. Since 77 is a function 
only of Aq/Aq [see Eq. (|A7p ] we can write all expressions 
for the phase diagram in terms of Aq instead of the actual 
transverse field H^. The transition line as a function of 
Aq is shown in Fig. Ofor ^0 ^ ^ (ic, for x^ 1 in the 
LiHoxYi_xF4 system); this diagram simply depicts the 
relation 1121 



Te(Ao) - 77(Ao)'r,(0) , 



(20) 



A. Classical Ising limit 

As shown above, if H±^ is small, the transverse hyper- 
fine interactions play a minor role (the A„j arc small), 
and the only effect of the longitudinal hyperfine inter- 
actions is to give a rather strong renormalization of the 
longitudinal dipolar interaction between the Ising dou- 
blet spins rj. The problem in this Ising limit (neglect- 
ing the transverse terms) was studied previously [331, but 
only for x= 1. We give a treatment here for all dilutions, 
and we also assume that Aq is arbitrary - surprisingly, 
the hyperfine interactions cannot be neglected even when 
^0 < Vq. 



If we now define e = (Tc — T)/Tc one finds that for 
Ao/^o,e < 1 (ie., smaU H±_ and T ^ Tc) the phase 
transition line Ac(r) obeys the relation 

A, =Aov^. (21) 

At T = Tc/2 one finds that A^ w Aq. When Aq/Ao > 1 
there is still a finite remnant polarization of the spin, 
and an ordered state at T = 0; the transition line obeys 
the relation Ac ~ Aq^J (Vq/T). This is quite different 
from the TFQI model (jT]), where for Ao > Vq the system 
becomes a paramagnet, and a T = quantum critical 
point is observed. 
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2. renormalized Ising model for arbitrary Ao /Vo 

We now relax the condition Vq <C so that all hyper- 
fine levels have to be included (however we still assume 
nuclear spin flips are blocked). 

For this case one can treat the Hamiltonian (jAl[) us- 
ing mean field theory; it then reduces to the mean field 
effective Hamiltonian 

i/MF = 5](wr-ffrK-E^"^r (22) 

i i 

where the site-dependent mean field is 

H:=T.^^(^D (23) 

3 

and h = luq. 

Since (|Aip is equivalent to the classical Ising Hamilto- 
nian (|A13p , the mean field Hamiltonian is equivalent 
to the mean field version of (|A13p . given by 

^ ^ '^im {^un ~t~ ^ini^im^ (24) 
irn 

I 



where now the mean field is 



Eim — y ] njm' Vim, im' i^jm' ) ■ (25) 
jm' 



The mean field theory in the form was solved some 
time ago [3^ for the homogeneous case (where (rj) is in- 
dependent of J, ie., the mean field is the same at all 
sites), and applied to the ferromagnetic LiHoF4 system 
(ie., when x= 1). 

In this section we extend this mean field approach to 
cover all values of x, including the spin-glass regime, 
by allowing the local mean field to vary from site to 
site. In order to allow easy comparison with the pre- 
vious work^^, we do this starting from the Hamiltonian 
in the form rather than ([M)) . An explicit derivation, 
given in App. |B1 results in the self consistent equation 



E 



m (/i2m2+Ag)3/2 



sinh 



cosh [P^jK^m? 



Ag) 



(26) 



When Aq ^ Vq the solution of this equation reproduces 
the results of Sec. IIV Al and in particular the phase dia- 
gram in Fig. [5] Let us now consider the regime Vq ^ Aq. 
Expanding Eq. (|26p in small Ag one obtains the behavior 
of the transition line Ac(r) near Tj,(0): 

Ac{e)^VoV~e, (27) 

where e = (Tc — T)/Tc- Because the dipolar interaction 
now dominates, the results of the TFQI model are repro- 
duced near T^, and at Tc/2 we find Ac ~ Vq. Surprisingly 
however, the hyperfine interaction, although small, dic- 
tates the physics at low temperatures. Below a crossover 
temperature T* = Aq/Vq one reaches a regime where 
Ac > Vb, and then Ac(T) is given by 

A,{T)=Ao^iVo/T) {T<T*) (28) 

corresponding to a transition temperature Tc{Hj_) = 
VqAq/Aq. These are the exact same formulae found 
above for the case Aq ^ Vq. Thus, when Aq > Vb the 
system gains more energy from fluctuations than it does 
from the interaction. However, since H± cannot flip nu- 
clear spins, a small remnant magnetization proportional 
to Aq/ Aq allows ordering at low temperatures. 

In this mean held theory one thus finds 2 regimes. The 
first, when Vq ^ Aq and Aq ^ Vb, is the standard Ising 



1 




FIG. 6; The magnetization as function of transverse effective 
field Ao is plotted for the Hamiltonian (|22p with h = 2Ao/7 = 
1, Vo — 20. Note that for Ao > Vb a remnant magnetization 
of magnitude h/ Aq is present. 



picture: at Aq = the spins are in either state |t) or 
and the electronic degrees of freedom order. For fi- 
nite Aq ^ Vb the spins fluctuate to the excited state 
at energy Vq. However, when Aq 3> Vq and/or in the 
whole parameter regime for Aq ^ Vq the physical pic- 
ture is different: the relevant single Ho Ising states are 
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An 




T 



(a) 



can be corrected for T > O.IK by choosing the scahng 
parameter better. However, in order to obtain a good 
fit with the experimental phase diagram at the lowest 
temperatures one has to take into account the transverse 
hyperfine terms [ij. For x< 1 the off-diagonal dipolar in- 
teractions have to be included as well. These interactions 
are considered next. 



B. Effect of transverse hyperfine interaction 



80 
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(b) 

FIG. 7: (a) The phase diagram of the Hamiltonian p2|) for the 
regime < Vo, with h = 2Ao/7 = l,Vo = 16. (b) Focusing 
on low temperatures, we compare the phase line in (a) (dashed 
line) to a system with Aq ^ Vq. We take ft = 4, Vb = 1 to 
have the same value for Ao\/\^, so that the low temperature 
divergence of the critical field is the same for the two systems. 



3|) , and the phase 
dependent infer- 



tile electro-nuclear states ||)-), Eq. 
transition line is dictated by their H± 
action, as discussed in Sec. lA II 

These two physical pictures arc best illustrated by the 
value of the magnetization at T = 0. For Aq ^ Vq, 
(X ?7, given in Ea. (|M)) . For Vq > and Ao < 
Vo, expanding Ea. (|B4[) in Aq/Vq, one sees that = 
1 — A§/(2Vo^), showing that the excitation energy is Vq- 
However, when Aq > Vq, the hyperfine energy dictates 
the magnetization, which is given by Mz « A/Vq for 
Ao = Vn and Ar = Aq/Aq for Aq > Vq (see Fig. El- 

In Fig[7] we plot the phase diagram of the mean field 
Hamiltonian ([221) as a function of T and H± for Vq ^ Aq. 
In the low T regime one can compare this with the phase 
diagram of a system with Aq ^ Vq and a similar value 
of Ao^/Vo. As expected from ([28l), for T <^ T* the two 
systems have the same behavior. 

In Ref. [s^l a similar phase diagram was calculated for 
LiHoF4 and compared to experiment [l| . This compari- 
son was made by rescaling the theoretical curve to agree 
with the experiments at the lowest temperature. How- 
ever the condition H± <C / is then not well satisfied 
at criticality, and the transverse hyperfine interactions 
are important. By forcing the theory and experiment to 
coincide in the regime where the theory is not applicable, 
a discrepancy with experiment over the whole tempera- 
ture range is obtained (sec Fig. 1(b) in Ref. [33| ) . This 



Independent of the ratio Aq/Vq, for Aq ^ Aq,Vq we 
found for ^^^yp = 0, which is equivalent to fJo ^ oo, that 
Ac = ^OA/(Vb/T), diverging as T ^ 0. This pathology 
arises because we need to include the transverse hyperfine 
terms. With H^y^ ^ the splitting A becomes apprecia- 
ble for ij,bH± « Aq « rio (in this regime ij,bH± « Aq), 
while Veff = VqAI/AI < Aq. Thus, for Aq^Vq ^ ^q a 
quantum phase transition is obtained within the regime 
of the applicability of the Hamiltonian (|A23[) . The di- 
vergence of Ac is rounded, as we schematically draw in 
solid line in Fig. [5l Similar rounding off of Ac occurs for 

Vq » Aq, Fig. m 

As mentioned above, the condition Aq <C i^o is not that 
well satisfied in the LiHoxYi_xF4 system. Still, at low x, 
where Aq ^ Vq, the phase transition occurs within the 
regime of applicability of the Hamiltonian (|A23[) : recall 
that in FiglHwe plot Vefj oc {Jz)^ for x= 0.167, taking 
Veff{H± = 0) to equal the value of Tc = 0.13K. The 
value of H± where A « Veff is smaller than H'^. This 
is true for all smaller dilutions x as well. At x= 0.167 
one expects the quantum phase transition to occur at 
H± w 2T, where A « Veff- Thus, three energy scales 
govern the phase transition. The spin-spin interaction 
Vq dictates Tc at zero field, the hyperfine interaction 
dictates the phase diagram at finite H±, and the larger 
anisotropy scale fig dictates the position of the quan- 
tum critical point, since quantum fluctuations only be- 
come important when (tiA'B-ff-Ll « f2o- It is for this 
reason 1 12l| that in LiHoxYi_xF4 it is much easier to dis- 
order the ordered phase thermally, rather than quantum 
mechanically [1], specially when x^ 1. 

For a; ^ 1 one can calculate the phase diagram for 
LiHoxYi_xF4 including all hyperfine terms using the 
effective Hamiltonian derived in Sec. IIIII However as 
we have seen, at higher temperatures this Hamiltonian 
breaks down at quite low transverse fields, because of the 
mixing of higher levels (cf. Fig. [4]). Therefore, for larger 
X, where the dipolar interactions arc stronger, quantum 
criticality occurs at H± > H*^ , where all the nuclear lev- 
els are well mixed. This is the case for x= 1 where a 
quantum phase transition is observed at 4.9T[i[. 

We therefore adopt a different approach, which covers 
all values of x, and calculate the phase diagram numer- 
ically, including both ^^fyp and i?hyp, starting from the 
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2.5r 




in which the single spin Hamiltonian is exact; and we 
then treat the interactions in mean field approximation, 
i.e., we assume 
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FIG. 8; (a) The phase diagram for x= 0.167 (thick lines) and 
x= 0.045 (thin lines) as a function of H± and T. The solid 
(dashed) lines are calculated with (without) the offdiagonal 
dipolar interactions. The low T, H± = phase is believed to 
be a spin glass for both dilutions, with a crossover at finite 
H± between a quasi spin- glass and a paramagnet. Note the 
"re-entrant" behavior predicted when the off-diagonal dipo- 
lar interactions are included (as seen by the crossing of the 
phase lines at H± ~ 1.65 T); see text, (b) The same phase 
diagram as in (a) now plotted as a function of Ao,T, to al- 
low comparison with experiment]^ at x= 0.167. Filled and 
empty circles, taken from Fig.l of Ref. [^, denote the PM-SG 
crossover from dynamical measurements and nonlinear sus- 
ceptibility respectively. With the inclusion of the off-diagonal 
dipolar interactions, good quantitative agreement is obtained. 



Hamiltonian 



gjHBH±Jj 



E 



TJZZ TZ TZ 



UMF{Jj)Jj 



(30) 



(see also App. |B|) . One of the central results of this paper 
is that the single atom Hamiltonian dictates much of the 
physics of LiHoxYi_xF4, and we shall see that the mean 
field approximation to the interactions has only a small 
effect on the results. 

The phase diagram for x= 0.167 and x= 0.045 is drawn 
in dashed lines in Fig [51 Comparing this calculation 
for x= 0.167 with experiment, we see it naturally ex- 
plains why it is much harder to disorder the spin-glass 
phase quantum mechanically than thermally. Going to 
x=0.045, we see the reduction in Tc is (x x, while the 
reduction in is much smaller, as can be anticipated 
from the requirement Vefj ~ Aq (Fig. S]). 

However, the agreement with experiment is still not 
perfect for small x. For x= 0.167 one obtains a larger 
critical field at T = and a qualitatively different behav- 
ior near Tc(0). As was discussed in Ref. [l^ . these dif- 
ferences can not be attributed to the mean field approxi- 
mation, but testify to the inadequacy of the Hamiltonian 
This is since the behavior near Tc(0) should fol- 
low Eq.(ini)(for x=: 0.167 the condition Vq < Aq is well 
satisfied), and the values for Ac obtained at the lower 
temperatures in the experiment Q necessitate the exis- 
tence of appreciable quantum fluctuations at H± fa IT, 
which contradicts the results shown in Fig[4l To explain 
things we now finally turn to the non-diagonal dipolar 
terms. 



C. Random non-diagonal dipolar terms 

To account for the experimental phase diagram, one 
has to include the dependence of the effective field on the 
offdiagonal terms of the dipolar interaction[13|. These 
add an effective random longitudinal field, and in the 
spin-glass regime also enhance the effective transverse 
magnetic field, as is explained in App. \^ Sec. lA 31 The 
random longitudinal field is crucial in dictating the na- 
ture of the phase at finite H±, as it destroys long-range 
SG order M, [Tgj : however, at least for x^ 1 it does not 
strongly affect the position of the phase line, because (i) 
the effective random longitudinal field is zero at H± — 0, 
and is small for H± <^ flo/jiB', (ii) it is random in sign, 
with only a small effect on the typical interaction; and 
(iii) at large H±, where 7f is appreciable, the crossover 
to the paramagnetic phase depends only weakly on Vq, 
as can be inferred from FigHl 

Thus, in calculating the phase diagram we neglect the 
random longitudinal fields, and consider only the en- 
hancement of the effective transverse field by the off- 
diagonal dipolar interactions. This enhancement depends 
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on X and H±; here wc follow Rcf.[l3| in neglecting the de- 
pendence on H±. We further assume that this enhance- 
ment is proportional to x when x^ 1, i.e. we write a 
total transverse mean field H± = H± + Hj^ , with H'l oc 
x. Note that this mean field is just the average of the 
transverse field H,]^ that we discussed in section Sec. lA 31 . 

This leads to a satisfying quantitative agreement with 
the experimental phase diagram at x= 0.167 (see Fig|S|). 
Tc(0) (X Vb (X x, while HI at the T = transition de- 
pends mainly on the energy scale JIq- For x< 0.167, HI 
should therefore change only slightly with dilution. The 
dilution dependence of H^ is a result of two effects. First, 
since the transition occurs when A « Vqj there is a slow 
decrease of HI with x (slow because A varies rapidly 
with H±; see FigU]). Second, Hf^ ex x, and giving a fur- 
ther reduction oc x in H^. This leads to the interesting 
prediction that for low enough x there will be an increase 
of HI with decreasing x, so that H'^ has a minimum at 
some x; this is seen in our figure by the crossing of the 
phase lines (see Fig|8]). In analogy with the re-entrant 
behavior one sees in some systems on variation of an ex- 
ternal field, we can call this a prediction of a kind of 
're-entrance' as a function of concentration x. 

It is interesting that the combined effect of the hyper- 
fine interactions and the transverse dipolar interactions 
leads to this re-entrant behavior. Even though the effect 
of the transverse dipolar interactions is only a weak ef- 
fect compared to that of the hyperfine terms, it is just 
enough to tip the system into re-entrance. Note however 
that without the much stronger hyperfine effect on the 
phase diagram, this would not have happened. We re- 
mark again that we do not think that it is possible to ex- 
plain the phase diagram without incorporating the hyper- 
fine terms (eg., by including only dipolar interactions [l8|; 
cf our discussion in section HIT)) . 



D. Nature of the low temperature phase 

As we have seen it is possible to derive an accurate 
phase diagram without saying too much about the nature 
of the phases themselves. In fact the nature of the low- 
T phases of LiHoxYi_xF4 has been rather controversial 
in recent years. Here we would like to outline several 
rather important implications of our results. We divide 
our discussion between the zero transverse field case and 
the case of finite H± . 



1. Zero Transverse field 

At all dilutions, the LiHoxYi_xF4 system is para- 
magnetic at high temperatures. However, as men- 
tioned above, at low temperatures the phase of the sys- 
tem is dilution dependent. It is well established both 
experimentally and theoreticallv[33j that for x> xp 
the system orders ferromagnetically at low temperatures, 



where values for xj? are in the range 0.2 — 0.5. However, 
at low dilutions the nature of the phase is controversial. 
Theoretically, it is argued that a spin-glass phase should 
exist at all dilutions x^ 1 3J| . Experimentally, it was 
argued that at x= 0.167 the system has a low tempera- 
ture glass phaseQ while for x= 0.045 the experimentQ 
revealed a very intriguing and yet unexplained behavior 
of the imaginary part of the susceptibility, in which its 
width in the frequency domain narrows as temperature 
is lowered, and therefore received the name "anti spin- 
glass". Recently, however, these results were challenged 
by Jonsson et al., [3^1, who claim for x= 0.167,0.045 
that there is no phase transition to the spin-glass phase. 
Furthermore, their analysis suggests that the system at 
the above two dilutions exhibits similar characteristics. 
A similar controversy arose regarding the specific heat 
of the system and its consequences regarding the nature 
of the phase at x= 0.045[l3, [H. Note, that it is dif- 
ficult to reach equilibrium conditions both experimen- 
tally, near the transition [s^, [sj, and numerically, using 
Monte Carlo [13, and therefore further studies will be 
useful in resolving the low temperature phase of the di- 
luted LiHoxYi_xF4. 

Our analysis above does not depend on the precise na- 
ture of the ordered phase, and therefore can not lead to 
definite conclusions regarding this question. However, 
since the Hamiltonian ^ gives a comprehensive descrip- 
tion of the system down to a few mK, some clarifying 
statements based on our analysis can be made. 

(i) the only differences between the LiHoxYi_xF4 com- 
pounds at x= 0.167 and x=: 0.045 are the strength of the 
dipolar interaction, both in the magnitude of the typ- 
ical terms and in the distribution due to randomness. 
All the single Ho properties, which, as discussed above, 
dictate much of the physics, stay unchanged. Thus, we 
have every reason to believe that at x= 0.045 the equi- 
librium low temperature phase is also a spin-glass. How- 
ever, as is shown in Fig[51 its Tc(0) is reduced to roughly 
SbrnK , and according to this the experiments at this 
dilution!^, [TtI . [ssI . l36l | were done in the paramagnetic 
regime. 

(ii) as we show above, the dynamics of the system 
at low temperatures is significantly slowed down by the 
couphng to the nuclear spins (see also Ref.js^). In- 
deed, thepeculiar features in the spin susceptibility at 
x= 0.045[8j were obtained as temperature was reduced 
to below 150mK. At this temperature the higher nu- 
clear spin levels start to be depleted, and all but few 
of the Ho atoms are in either state |t, —7/2) or |J,, 7/2). 
Thus, the system can not take advantage of the much 
faster transitions between the higher nuclear spin states 
(see FigU]), and the dynamics slow down appreciably. 
The data of Quilliam et al.js^, showing that the peak 
in the specific heat occurs in a similar temperature for 
x= 0.02,0.045,0.08 supports the view that single spin 
physics, and in particular the hyperfine interactions, are 
significant in the interpretation of the experiments in 
these dilutions. 
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(iii) In Ref. [l7[ it was argued that for x= 0.045 the 
internal transverse field resulting from the offdiagonal 
terms of the dipolar interaction stabilize a low tempera- 
ture spin liquid state. It was further argued there that 
this is correct also for transverse fields which are effec- 
tively reduced by a factor of 10"*. The analysis of Ref. [l3| 
was done in the electronic degrees of freedom. However, 
in the regime relevant to the experiment Q the effective 
Hamiltonian (|A23p is valid, with zero random longitudi- 
nal field. Therefore, the analysis should be done consid- 
ering the electronuclear degrees of freedom, within the 
framework of the Hamiltonian (|A23[) . In particular, the 
effective transverse field due to the offdiagonal dipolar 
interactions at H± = is much smaller than the values 
considered in Ref. [l3|, as can be inferred from the log 
scale graph in Figdl 



2. Finite Transverse Field 

Turning now to non-zero H±, we note first that there 
is a crucial difference between the FM phase and the SG 
phase, if the latter exists at H± = 0. In the FM regime, 
the lower critical dimension dc — 2 (cf. ref. [Ill), and in 
3D the FM phase is stable to a small random field. Thus 
we expect that the FM phase, which exists for large x, 
will survive at finite H±. 

However, if one supposes that for intermediate x one 
has a SG phase at H \ = 0, then the critical dimension is 
dc = oo (cf. refs.[3^,|43|), and so the long-range SG order 
should be destroyed by an infinitesimal random field [H, 
[3. As is well-known, this means that the system will 
no longer be a homogeneous SG, but instead domains of 
finite size will be created - each one will have internal SG 
order but the order will be uncorrelated between different 
domains. The correlation length ^, which is essentially 
the domain size, is given bvfTa. [l^ 



/ no 



(3/2)-e^ 



(31) 



where 6^ ~ 0.2 is the stiffness exponent [391 140|. Essen- 
tially the system is able to gain energy from the random 
field by creating domains. Referring to Eq. (|A27[ ). we see 
that this energy gain is a result of the two terms in the 
numerator contributing with the same sign, i.e. an effec- 
tive enhancement of the transverse magnetic field. 



V. EXPERIMENTAL CONSEQUENCES 

In this paper we have derived new results regard- 
ing the single particle properties of the Ho ion in the 
LiHoxYi_xF4 systems, as well as the phase diagram. We 
have also addressed the regime where Aq <^ Vq, which 
is not applicable to the LiHoxYi_xF4 system but is the 
more abundant regime in general. In this section we ad- 
dress the relation between our results and possible exper- 
iments. 



A. Single spin properties 

A central result of this paper is the derivation of the 
low energy effective Hamiltonian for the LiHoxYi_xF4 
system, as a generalized Ising Hamiltonian in the elec- 
tronuclear degrees of freedom (fT5|) . This effective Hamil- 
tonian is completely determined by the effective random 
fields 7j^, and the single ion parameters Cm, Am, and 
r],n [see Eq. (|A16|) ]. the latter determining the effective 
spin and therefore the effective spin-spin interaction. Be- 
low we suggest magnetic resonance and fj.SK experiments 
that can measure the single ion parameters directly, and 
verify the mechanism leading to the enhancement of the 
effective transverse field and the emergence of an effec- 
tive random longitudinal field. With regard to magnetic 
resonance experiments, we give explicit quantitative pre- 
dictions for the Rabi frequency of excitations to various 
levels. We interpret these predictions in terms of the 
calculated entanglement entropy of the ground state as 
function of transverse field, in agreement with our anal- 
ysis in Sec. Ill CI 



1. Magnetic Resonance Experiments 

One obvious way of probing the low energy properties 
of the LiHoxYi_xF4 system is via magnetic resonance 
experiments. Specifically, such experiments can be used 
to quantify £„, Am, and the nature of the wave func- 
tions as function of Hj 
frequencies, given by 



In Fig. Ufa) we plot the Rabi 



26^ I (1 I g.jfiBJz + gNfJ-Nh I 2) 



(32) 



for a magnetic resonance transitions between the ground 
state and excited states as function of transverse field 
for an AC field along z. Only the lowest 16 levels are 
considered. From symmetry, only transitions to antisym- 
metric states (plotted dashed in Fig [31 see also discussion 
in Sec. Ill C[) are possible. For H± the only allowed 
transition is to the first excited state, in agreement with 
the form of the electronuclear states in Eq. (|A3[) . in terms 
of which the effective Hamiltonian ^ is written. Finite 
Rabi frequency to other states is allowed at 7^ 0, con- 
sequence of the mixing of the states ([A3p resulting from 



Jfj^^and is thus larger for larger |m|. As discussed in 
Scc lII CI with increasing the electronic and nuclear 
spins disentangle, and the electronic state at high fields is 
approximately the symmetric state (|t) + |i))/V2 for the 
lower 8 states and the antisymmetric state (|t) ^ U))/V2 
for the upper 8 states. For this reason, except for the first 
excited state at small H^, the Rabi frequency is larger to 
the states in the upper group. At large Hx The Rabi 
frequency to state 9 (the lowest level in the upper group) 
dominates, in agreement with the picture (see Sec. Ill C[) 
that levels in the lower and upper groups of 8 have similar 
nuclear states, respectively. 



13 



100 



80 



^ 60 



S 40 



20 







|m|=7/2 '^\^ 




\ / 




\ / 

y 


- 


/ \ 


^>^^4rTij=3/2 








^II;;::s«ii=5/2_ 



1 2 

(Tesia) 

(a) 




(TesIa) 

(b) 



FIG. 9: (a) Rabi frequencies of magnetic resonance transitions 
to electro- nuclear levels within the space of the electronic Ising 
doublet, for oscillating field along z with amplitude ImT, as 
function of a transverse magnetic field. All transitions are to 
anti-symmetric states (dashed blue levels in Fig|3}. Dashed 
(solid) lines correspond to levels within the lower (upper) 8 
states, (b) Entanglement entropy of the ground state as func- 
tion of _ffx (solid black line). Plotted for comparison are the 
normalized Rabi frequency to level 2 (first excited state, blue 
dot dash line) and the subtraction from unity of the normal- 
ized Rabi frequency to level 9 (green dashed line). 



The disentanglement of the electronic and nuclear 
states can be quantified by calculating, as function of 
Hx^ the entanglement entropy ~Tr{pj logp^), where 



re-plot the Rabi frequency to the first excited state scaled 
to 1 at 77a; = 0. In dot-dash line we plot 1 — t'RabiCS), 
where !>Rabi(9) is the scaled Rabi frequency to level 9. Al- 
though not exact, we see that the diminishing of z^Rabi(2) 
and the emergence of i^Rabi(9) with increasing field is a 
measure of the (dis)cntanglcmcnt of the electronic and 
nuclear spins. A naive conclusion from the above would 
be that for large H±^ the nuclear spins decouple from 
the electronic spins, and therefore the effective Hamilto- 
nian ([1]) is recovered. The fact that for x= 1 the soft 
mode is gapp ed near the quantum phase transition (at 
4.9T)(l^ Il5| suggests that this simplified model is not 
suitable also in this regime. 

In Figlin] we plot the Rabi frequency as function of 
for an AC field in the x direction. The relevant ma- 
trix element is then (1 | gjfJ-BJx + gNfJ-Nlx \ 2). The 
operator changes the z component of the electronic 
spin. The relevant matrix element is then proportional 
to the amplitude of | Fj) in the ground state. This ampli- 
tude, resulting from the transverse hypcrfinc interaction, 
is small, ~ O(^o/^^o) ~ 10~^, and at = is finite 
only for the state with | \— 5/2. This is why the in- 
tensity for a longitudinal AC field is so much larger than 
that for a transverse AC field (compare Fig. EJa) and 
Fig. [Tn]) . For an AC field along z, only transitions to 
symmetric states (plotted sohd in Fig[3]) are possible. 

From the picture of the system without the transverse 
hyperfine interaction, where the 16 states at zero field are 
eigenstates of , one might expect that for a transverse 
AC field the dominant intensity would come from the 
nuclear operator. Surprisingly, it is the electronic oper- 
ator that dominates the magnetic resonance experiment. 
This is because Aq/^Iq 3> ^n/mb- Thus, although the 
levels are predominantly nuclear spin levels, the relevant 
experiment is basically an ESR experiment. 

All the results above are valid for single Ho ions, and 
can be checked in very dilute samples, where interactions 
are negligible. For larger x, in the SG regime, the inter- 
play of the offdiagonal dipolar interactions and the ap- 
plied transverse field results in an effective enhancement 
of the transverse field and the emergence of an effective 
random field[l3, (see details in Ap p. \M Sec. [X3|) . 
This result, shown in Sees . IIV Ul IIV D 2] to be crucial for 
the structure of the phase diagram in the SG regime 
as well as for the nature of the phase itself[l6l.[l9t. can ac- 
tually be verified by measuring e.g. i^Rabi(2) as function 
of H± for different x. The offdiagonal dipolar interac- 
tions should lead not only to a dispersion in A7/2J but 
also to an x dependent shift upward of its mean value. 
This shift can be checked against our approximation in 
Sec. UVCl 



p,=Y,{M\g.s.){g.s.\M) (33) 

2. i-iSR Experiments 

is the reduced density matrix in the subsystem of the nu- 
clear spin. The entanglement entropy is shown for the According to Eq. (jA6[) . the single Ho spin moment at 
ground state in solid line in Fig.[5I^b). In dashed line wc low temperatures oc 77, and so decreases with H±. Such 
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FIG. 10: Rabi frequency of magnetic resonance transitions 
for an AC field in the x direction. Details are as in Fig[9] 
only here the transitions are to symmetric states (solid green 
levels in FiglS]). 



a field dependence of the individual magnetic moments 
could be directly measured using //SR. The magnetic field 
at the muon site is proportional to the magnetic moment 
size of the material, which, in a diluted sample, is given 
by the nearest Ho ion. Such a measurement should be 
done at dilution x<C 1, both because our prediction is 
for the regime where > Vb, and because then the 
contribution from more distant Ho ions will be smaller. 



B. LiHoxYi_xF4 Phase diagram 

In Scc lIV CI We have given explicit quantitative predic- 
tions for the phase diagram of the LiHoxYi_xF4 system 
as function of x. Our results for x=0.0167 are in good 
agreement with experiment, and our predictions regard- 
ing the phase diagram at x= 0.045, and the re-entrant 
cross-over field as function of dilution, can be checked 
experimentally in a straight forward way. However, here 
we would like to suggest an experiment that would di- 
rectly probe the significance of the hyperfine interactions 
in dictating the phase diagram of the LiHoxYi_xF4 sys- 
tem at low X. 

We use the above result, predicting that at ~ 
2 — 3T, the magnetic resonance intensities for transitions 
to levels with | |= 1/2,3/2 arc appreciable. Thus, one 
could in principle, by populating these states, change the 
critical field at low-T: a non-equilibrium occupation of 
these excited electro-nuclear levels would lead to stronger 
quantum fluctuations (cf. Fig. 3]), and therefore to a 
lower critical field. This opens up the rather fascinat- 
ing possibility of controlling the quasi-equilibrium phase 
diagram of the system by driving a steady-state non- 
equilibrium nuclear spin population. 



C. Limitations of the LiHoxYi_xF4 system 

The LiHoxYi_xF4 compound is a particularly useful 
test system: It is a well defined Ising system, with a 
doubly degenerate ground state, quantum fluctuation are 
easily tunable at moderate transverse fields, and x can be 
varied over a huge range. However, there are at least two 
limitations on this system, viz.: 

(i) The allowed values of the dipolar spin-spin interac- 
tion Vo, hyperfine interaction Aq^ and crystal anisotropy 
energy JIq, do not test the whole parameter range. Thus, 
eg., to observe our prediction that the hyperfine coupling 
dictates a diverging H^, at low T for cither Aq <C Vq or 
^0 ^ given that Aq, ^ ^^o (see Figs l5l7p . we need 
a system where the latter condition is well satisfied. 

(ii) For iJ_L = 0, the ground state is degenerate. Induc- 
ing quantum fluctuations coupling the two ground states 
requires a transverse field. However, the application of 
ff^ results in an emerging random field. As a result, 
the quantum phase transition between the SG and FM 
phases can not be seen as function of but only as 
function of a parameter that docs not break time re- 
versal symmetry, e.g. pressure jlGj. To observe such a 
transition one would need a system where quantum fluc- 
tuations between the Ising ground states are appreciable 
at iJ^ = 0. One would then have to tune the dilution 
so that at ambient pressure the typical spin spin inter- 
actions is of the order of the quantum fluctuations, and 
look for the transition as function of pressure. 



VI. CONCLUSIONS 

In this paper we have considered anisotropic quantum 
magnetic systems in which both dipolar and hyperflne 
interactions play a role. We have shown that the trans- 
verse field Ising model is not sufficient to describe such 
systems; instead, we have given a theoretical treatment of 
an Elcctronuclcar Quantum Ising model which can do the 
job. The hyperfine interactions set the scale for the field 
at the quantum critical point, even in systems in which 
the hyperfine interaction is weaker than the dipolar spin- 
spin interaction. We have given a detailed treatment of 
the LiHoxYi_xF4 compound, calculating the phase di- 
agram for all dilutions x, and giving explicit numerical 
results for x= 0.045 and x= 0.167. We explain the ex- 
perimental result that thermal fluctuations more easily 
destabilize the ordered phase than quantum mechanical 
fluctuations. Off-diagonal dipolar interaction terms are 
shown to reduce the transverse critical field ^ and a 
prediction for a non-monotonic critical field as a function 
of X is given. Experimental consequences of our results 
as well as possible measurements of the parameters of the 
effective Hamiltonian are discussed. 

We note that our results have wider implications in 2 
ways, which will be explored elsewhere. First, as just 
noted, they can be applied to many other dipolar quan- 
tum magnets. Second, the nuclear spins will clearly have 
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an even more profound effect on the dynamical proper- 
ties of these systems than on the phase diagram - indeed, 
the big surprise is quite how important they are for the 
thermodynamics, even when ^ 

Note added in proof - Our prediction for the rccntrance 
of the crossover field as function of dilution was coinci- 
dentally and independently discovered experimentally by 
Ancona- Torres et. al.[4l|. 

It is a pleasure to thank G. Aeppli, B. Barbara, B. 
Malkin. A. Morello, and J. Rodriguez, for useful discus- 
sions. This work was supported by NSERC in Canada, 
and by PITP. 



APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN 

In this Appendix we derive the effective Hamiltonian 
in Eq. (HT 



1. Ising-like terms 

We are interested in the 21 + \ lowest levels, shown 
for the 16 levels of the Ho ion in LiIl0xYi_xF4 in FiglS] 
The term in Eq.llSl), has zero matrix element be- 
tween any of these 16 low energy states, and thus must 
mix higher crystal field states. This results in contribu- 
tions of order Vb/iloi ^j/^^o ^ 1- We therefore begin by 
discussing in this subsection the longitudinal interaction 
terms, i.e. the Hamiltonian Hi = Hq + Hf^. 

In the subspace of the lowest 21+1 clcctronuclcar 
states, we have seen that Hi reduces to 



where | a), | b),\ a),\ h) are defined in Sec. Ill Bl and 



Cl 



C2 



(A4) 
(A5) 



Thus, as noted just above, the longitudinal part of the hy- 
pcrfinc interaction blocks quantum fluctuations between 
the relevant Ising states. 

It then follows that a transverse field H± <C ilo can 
only renormalizc the effective spin of what is just a clas- 
sical Ising system; one finds 



(4(i^±)> 
with (rl) = 



'/(4(0)) 



V^{cl-cl), (A6) 



Note that 



A 



^1 



1 - ^ (Ao « A,) 
7] = Ao/Ao (AoMo > 1) 



(A7) 



Absorbing this rcnormalization into the dipolar interac- 
tion, Hic (|A2[) reduces to: 



hL = 



E 



[H^ « Oo/mb) (A8) 
(A9) 



where Sj is a spin-half matrix operating on the states 
III-) and I-IJ-) of the j-th spin, such that s| jfl") =111), and 
s||4) = - 14), etc. 

Thus we have shown the equivalence of the two Hamil- 
tonians H\f, in \K2\ and ijjjj in (|A9[) . and both arc ap- 
plicable in the low-T limit kT <^ ujq, when H± <C i^o- 

For higher temperatures, luq < T <^ ^Iq, the Hamilto- 
nian (|A1[) has to be considered. Following the arguments 
above, only levels with the same mix (see Fig.[2|), and 
the generalization of Eqs. (jA3| )- (|A5p results in 



Let us first discuss this Hamiltonian for T Uq. Then, 
one can simplify the model to include only the levels 
with = ±1 (ic. with 4 = ±7/2 for LiHoxYi_xF4), 
described in the previous section. The effective Hamilto- 
nian then reduces to [131: 



ffle = 2AoE' 



-EA^ 



(A2) 



where acts upon the two electronic states IT),!!) at 
site i, Ui acts upon the two nuclear spin states with Iz = 
±7/2, and the coupling Aq — Iluq ~ 0.7 K. 

For H± the two low energy Ising doublet states 
are given bv[l^ 



1^) = Cl I a) + C2 I 6) 
m = Cl I a) + C2 I 6) 



(A3) 



|-f|-,m) = ci,„ |t,m) +C2m |i,m) 
H, -m) = cim ||,-m) +C2m |T,-"^) , (AlO) 



with coefficients 



Clm — QfmAo; 

C2m = Ctm ["ICJQ + \/ m'^i^o + Ag] (All) 



and 



A^ + {mujo + \lm?ujl + Ag)^ 



-1/2 



(A12) 



If we choose m = — /, these equations revert to (|A3[) - (|A5[) 
for the 2 lowest levels. 

We now define a set of pseudospin-1/2 degrees of free- 
dom {sim}, operating in the subspace spanned by the 
pair of degenerate levels | ff,™), | -IJ-, — m). In terms of 



16 



these pseudospins we then obtain, for the Hamihonian 
Hi defined above [Eq. (|Aip ]. the renormahzed effective 
Hamiltonian 

jm 

~ ''^im'>T'jm'Vijn,jm'i-^±)^im^jrn'i-^^'^) 

where we have defined a psedospin occupation number 
rijm such that ~ 1 , and defined energies 

em(ff±) = sgn{m)^{m^ujl + Al), (A14) 

and a renormahsation factor 

Vrn = (ci„ - cl,J. (A16) 

Thus again we obtain a generahzed classical Ising model, 
now involving the entire set {sim} of pseudospins. There 
are 2/ + 1 pseudospins per site (only one of which is 
occupied at any time), interacting with pseudospins at 
the other sites. 

We emphasize that the Hamiltonians (|A1[) and (|A13[) 
are entirely equivalent. The great advantage of (jAl3|) 
(and its low-T simplification in ()A9P ) is that the physics 
is correctly displayed, that of a classical Ising system; 
and this is done using the physically meaningful energy 
scales for this regime. 

2. Transverse Hyperfine terms 

Now suppose we neglect all dipolar interactions be- 
tween the electronic spins, but we now switch on the 
full hyperfine coupling, including the transverse hyper- 
fine term. The general effect of this is seen in a plot (Fig. 
[3]) of the 16 relevant eigenenergies for the Ho ion in the 
LiHoxYi_xF4 system, as a function of H±. 

As H± increases, levels separate into 2 groups of eight, 
given by symmetric and antisymmetric combinations of 
states with the same m. Each pair of levels IfT, m), |J| 
, —m) which is related by time reversal symmetry is then 
split by the combination of H± and H^^^. 

If we wish to write an effective Hamiltonian for this 
system in the original basis of 2J -I- 1 levels, we get a 
rather interesting result. After truncating the full Hd 
down to the 2 lowest electronic levels, we can write for a 
single ion (ignoring now all interactions between ions): 

j i 

where the matrix element fjj* = (./le*^'"'"^' |i) is defined 
between an initial state |j) = | T j;Xi(Ij)) before the 



electronic spin flips, and a flnal state |/) = | i j; x/(Ij")) 
after it flips, involving some initial and final nuclear spin 
wave- functions Xi(Ii) ^-^id x/(Ij)- The operator t/'^' is 
defined as 

j)_g.A[i,] ^ g^p_ dtH^^^{lj,t) (A18) 

and in general the (2/ -I- 1) x {21 + l)-matrix A[I], which 
operates in the Hilbert space of Ij, causes transitions 
between different hyperfine levels when the central spin 
(here the Ho spin) flips. A calculation of A[I] is actually 
quite lengthy, since typically it involves multiple transi- 
tions between the nuclear spin states. Thus the Hamil- 
tonian (jA17[) is in general a rather complicated object, 
taking the form of a (2/ -I- 1) x (2/ -I- l)-matrix acting on 
the lowest states of the system. For LiH0xYi_xF4 this 
means a 16 x 16 matrix, which is quite unwieldy. 

In the present paper we are only interested in the ther- 
modynamic properties of our system. This allows a con- 
siderable simplification, which we can see most simply 
by rewriting the Hamiltonian (jA17p for an isolated ion 
(ie., again ignoring dipolar interactions) in terms of our 
pseudospins, in the following form: 

H'J^ = Y.^°ni.S^r.^)+ ^ '^^W (^zm , W ) (A19) 

where the 'diagonal' terms have the matrix form 




(A20) 



in the basis where is diagonal, and the non-diagonal 
terms couple different pseudospins. 

We look first at the diagonal terms. The diagonal 
energies e„i (|A14p are just the eigenvalues of in 
(fTSll : in zero field one has from (fT3|) that e„i ~ mujQ 
for LiHoxYi_xF4, with splitting ^ 205 mK between ad- 
jacent pairs of levels in the Ho ion. At low H± the 
transition terms Am split each pair of degenerate m- 
states into the symmetric and antisymmetric combina- 
tions |±,to) = (I -f|-,m) ± I JJ., — to)/V2. The A™ are just 
the quantum fluctuation amplitudes between the eigen- 
states I 11, m) and | J|, — m) of the classical Ising system, 
induced by the transverse hyperflne coupling. 

Now consider the non-diagonal term 5Hmm'{sirm Sim') 
in (jAigp . The crucial point here is that when H^/VLq <C 
1. this term will be unimportant for the phase diagram 
because a non-diagonal couplings between different pseu- 
dospins Sim, Sim' involves a sequence of |to — m'| nu- 
clear flips. If we call these non-diagonal matrix elements 
Amm', then for H±/flo <C 1 Amm' ^ 'j-'o and so it can 
hardly affect the level spacing or any other thermody- 
namic properties. Thus, as far as the thermodynamics is 
concerned, we can get away with using the form H°j^ in 
dUni, when H±/no < 1. 

Summarizing, including the transverse hyperfine terms 
the effective low energy Hamiltonian for a single ion can 
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be written as 

Hf - - A,„(iJj_)sf,„] (A21) 

where both e,„ and depend on the transverse field. If 
we now include back the longitudinal dipolar interaction, 
we obtain a Hamiltonian 

-E.,™"™[e™ + A™(i?x,wo)<J (A22) 

which is the generalization of Eq. (jA13|) that now includes 
quantum fluctuations. 

At very low temperatures, T ^ luq (which, for x<C 1, 
include the whole phase diagram) only the 2 lowest 
electro- nuclear states are relevant, and the above Hamil- 
tonian reduces top^ 

^ = - E ^'(H^)s!s^ - A(iJx) ^ , (A23) 

which is just Eq. (jA9p with the addition of quantum fluc- 
tuations, and is the ENQI Hamiltonian given in Eq. 

We have still not quite finished; wc must finally add 
in transverse dipolar terms, which introduce one further 
modification to the effective Hamiltonian. 



-ffiong has a low-T ordered phase, which, depending on 
the dilution, is FM or a SG@,l,[3l. Let us denote either 
of the two symmetry broken ground states of the ordered 
phase by tpo- Introduciiig if trans as a perturbation lowers 
the energy of tpo by[l^ Elj H^l 

^(2) ii'o U^JN! + gjf^BH^ Jr?\ V'o) 

no ■ 

(A27) 

Thus, the offdiagonal dipolar terms add to the applied 
transverse field a term 



77W(r,) = ^ '\ (A28) 

This additional field is random in sign, and can enhance 
or decrease the quantum fluctuations generated by the 
applied field H± . Since enhancing quantum fluctuations 
reduces the energy of the system, configurations in which 

(r) 

Hx is in the direction of the applied field are energet- 
ically favorable. The terms in Eq. (jA27[) proportional 
to H]^ and to are independent of the spin configura- 
tion of the system. However, the cross term depends on 
the location and orientation of the s pins , and results in 
an effective random longitudinal field|l6l. [Tgl. [20| j^{H±) 
given by 



3. Non-diagonal Dipolar Interactions 

The 'non-diagonal' dipolar terms couple the Ho 
spins in higher order in the small parameter U[j /Hq, and 
so they have typically been neglected when discussing 
anisotropic dipolar systems, including LiHoF4. However, 
in Refs. [H, [ij, [l^ shown that terms ~ Jj^i can 

be rather important. For x= 1 these terms cancel by sym- 
metry, but not for < x < I, where even when H± = 
they induce quantum fluctuations [l3l. For H± ^ these 
terms can enhance or reduce the quantum fluctuations 
induced by the applied field. To quantify this effect let 
us consider the regime H± <C and write the original 
Hamiltonian in Eq.lg]) or Eq.((Tl]) as 

H — i^long + -ff trans (A24) 

where 

ffiong = Hcf + iihyp -I- J7dip (A25) 

and 

iftrans = " 9jmH±r + ^ C/^W^ (A26) 
i ij 

Here we neglect ff^yp, and all terms other than Ufj' in 
^dip which do not contribute in lowest order of pertur- 
bation theory [H, [3. 



7^ y^Uji = Uo, (A29) 

3 

where J = (J^) « 5 is the single spin magnetic moment, 
Ci is a random number with (ci) = and Var{ci) = c^(x) 
is dilution dependent. For (1 — x) <g; 1 we have c(x) = 
c'(l — x)[2^, with c' w 1. Note that the effective ran- 
dom field exists both in the SG[ll,[lil and in the FMjl^l 
regimes. The randomness is a result of the quenched 
disorder, so the interplay between the applied transverse 
field and the off-diagonal dipolar terms converts spatial 
disorder to randomness in the effective longitudinal field. 
Unlike the quantum fluctuation amplitudes A,„ , this ran- 
dom field is independent of m to zeroth order in wq/^^o- 
It is clear from the definitions of the additional trans- 
verse field ni^'' {vi ) that its actual values and distribution 
depends on the particular configuration adopted by the 
{^im,} i'^ the phase of interest. Thus this final term in the 
effective Hamiltonian actually depends on what state the 
system is in. In general this extra field simply renormal- 
izes the total field acting on the system; we have a new 
transverse field 

Ht = H^+H'^;\r,). (A30) 

This means that all the parameters in the effective Hamil- 
tonian that formally depend on the transverse field must 
now depend on rather than H±. This inevitably in- 
troduces some randomness in these parameters, from one 
site to another. However, as is discussed in Sec. IIVDI 
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in the SG regime, the system forms finite size domains 
that maximize the energy gain from the random field . 
This is actuahy done by having a finite average value 
for Hi^\ri) in the direction oi H±, thus increasing quan- 
tum fluctuations, and results in an effective enhancement 
of the applied magnetic field. With the addition of the 
effective random field and the effective enhancement of 
the transverse field we obtain the effective Hamiltonian 
in Eq. (UHl). 



APPENDIX B: CLASSICAL ISING LIMIT IN 
MEAN FIELD 



and the average magnetization of spin i is given by 



exv\-mn - - Aorf]]} • (B3) 



Equations (|B2p and (|B3p are generalizations of Eqs.(6) 
and (10) of Ref. [12 to the case where Vi is site dependent. 
Note that Vi is positive by definition. Assume now that 
the paramagnet-spin-glass phase transition occurs when 
the mean field equation gives a finite magnetization for 
spins at sites with some typical Vi — Vq- Then, defining 
Mz = (r^) one obtains 



In this Appendix we derive Eq. (|^5|) . Consider a given 
site i, with local longitudinal field Hi. We assume that as 
we cross the transition line into the spin-glass phase the 
expectation value of each spin grows at the same rate. 



I.e. 



with Vi - 
comes Z 



It then follows that 

H^ = v(T,n 



(Bl) 



^jVij^aij. Now the partition function be- 
n,; Zi with 



Z. = Tr{exp[-f3[{hl! " HHr^ - Aorf]]} (B2) 



j:^cosh[l3H{m)] 



(B4) 



where H{m) ~ [hm + V^Mz)"^ + Aq is the total mag- 
nitude of the mean field. This equation allows us to ob- 
tain the phase diagram for any ratio of Aq/Vq, keeping 
H± ^ rto/fiB- Near the transition line (so <^ 1) 
and following Banerjee and Dattagupta[33|, we expand 
Eq.JHU in and obtain Eq. ([26l) . 



[1] D. Bitko, T. F. Rosenbaum, and G. Aeppli, Phys. Rev. 

Lett. 77, 940 (1996). 
[2] D. H. Reich, B. Ellman, J. Yang, T. F. Rosenbaum, 

G. Aeppli, and D. P. Belanger, Phys. Rev. B 42, 4631 
(1990). 

[3] T. F. Rosenbaum, J. Phys. Condens. Matter 8, 9759 
(1996). 

[4] W. Wu, B. Ellman, T. F. Rosenbaum, G. Aeppli, and D. 

H. Reich, Phys. Rev. Lett. 67, 2076 (1991). 

[5] W. Wu, D. Bitko, T. F. Rosenbaum, and G. Aeppli, Phys. 

Rev. Lett. 71, 1919 (1993). 
[6] J. Brooke, T. F. Rosenbaum, and G. Aeppli, Nature 413, 

610 (2001). 

[7] J. Brooke, D. Bitko, T. F. Rosenbaum, and G. Aeppli, 

Science 284, 779 (1999). 
[8] S. Ghosh, R. Parthasarathy, T. F. Rosenbaum, and G. 

Aeppli, Science 296, 2195 (2002). 
[9] R. Giraud, W. Wernsdorfer, A. M. Tkachuk, D. Mailly, 
and B. Barbara, Phys. Rev. Lett. 87, 057203 (2001). 
[10] R. Giraud, A. M. Tkachuk, and B. Barbara, Phys. Rev. 

Lett. 91, 257204 (2003). 
[11] B. Barbara et. ai, J. Magnetism and Magnetic Materials 

272-276, 1024-1029 (2004). 
[12] M. Schechter and P. C. E. Stamp, Phys. Rev. Lett. 95, 
267208 (2005). 

[13] N. V. Prokof'ev and P. C. E. Stamp, Rep. Prog. Phys. 
63, 669 (2000). 

[14] H. M. Ronnow, R. Parthasarathy, J. Jensen, G. Aeppli, 
T. F. Rosenbaum, and D. F. McMorrow, Science 308, 
389 (2005). 

[15] H. M. Ronnow, J. Jensen, R. Parthasarathy, G. Aeppli, 



T. F. Rosenbaum, D. F. McMorrow, and C. Kraemer, 
Phys. Rev. B 75, 054426 (2007). 
[16] M. Schechter and N. Laflorencie, Phys. Rev. Lett. 97, 
137204 (2006). 

[17] S. Ghosh, T. F. Rosenbaum, G. AeppU, and S. N. Cop- 
persmith, Nature 425, 48 (2003). 

[18] S.M.A. Tabei, M.J. P. Gingras, Y.J. Kao, P. Stasiak, and 
J.Y. Fortin, Phys. Rev. Lett. 97, 237203 (2006). 

[19] M. Schechter, P. C. E. Stamp, and N. Laflorencie, J. 
Phys.: Condens. Matter 19, 145218 (2007). 

[20] M. Schechter, Phys. Rev. B 77, 020401(R) (2008). 

[21] P. B. Chakraborty, P. Henelius, H. Kjonsberg, A. W. 
Sandvik, and S. M. Girvin, Phys. Rev. B 70, 144411 
(2004). 

[22] S. Bertaina, B. Barbara, R. Giraud, B. Z. Malkin, M. 

V. Vanuynin, A. I. Pominov, A. L. Stolov, and A. M. 

Tkachuk, Phys. Rev. B 74, 184421 (2006). 
[23] N.V. Prokof'ev, P.C. E. Stamp, J Low Temp. Phys. 104, 

143 (1996). 

[24] N. V. Prokof'ev and P. C. E. Stamp, Phys. Rev. Lett. 

80, 5794 (1998). 
[25] P.C.E. Stamp, I.S. Tupitsyn, Phys. Rev. B69, 014401 

(2004). 

[26] P. E. Hansen, T. Johansson, and R. Nevald, Phys. Rev. 
B 12, 5315 (1975). 

[27] G. Mennenga, L. J. dejongh, and W. J. Huiskamp, J. 
Magn. Magn. Mater. 44, 59 (1984). 

[28] G. S. Shakurov, M. V. Vanyunin, B. Z. Malkin, B. Bar- 
bara, R. Yu. Abdulsabirov, and S. L. Korableava, Ap- 
plied Magnetic Resonance 28, 251 (2005). 

[29] D. M. Silevitch, D. Bitko, J. Brooke, S. Ghosh, G. Aepph, 



19 



and T. F. Rosenbaum, Nature 448, 567 (2007). 
[30] S. M. A. Tabei, F. Vernay, and M. J. P. Gingras, Pliys. 

Rev. B 77, 014432 (2008). 
[31] K. Andres, Phys. Rev. B 7, 4295 (1973). 
[32] V. Banerjee and S. Dattagupta, Phys. Rev. B 64, 024427 

(2001). 

[33] A. Biltmo and P. Henelius, Phys. Rev. B 76 054423 
(2007). 

[34] M. J. Stephen and A. Aharony, J. Phys. C 14, 1665 
(1981). 

[35] P. E. Jonsson, R. Mathieu, W. Wernsdorfer, A. M. 
Tkachuk, and B. Barbara, Phys. Rev. Lett. 98, 256403 



(2007). 

[36] J. A. Quilham, C. G. A. Mugford, A. Gomez, S. W. Ky- 
cia, and J. B. Kycia, Phys. Rev. Lett. 98, 037203 (2007). 

[37] B. Barbara, Phys. Rev. Lett. 99, 177201 (2007). 

[38] Y Imry and S. K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 

[39] D. S. Fisher and D. A. Huse, Phys. Rev. Lett. 56, 1601 
(1986). 

[40] D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 386 (1988). 
[41] C. Ancona-Torres, D. M. Silevitch, G. Aepph, and T. F. 
Rosenbaum, larXiv:0801.2T8n 



